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Abstract 

The  laminar  convective  heat  transfer  and  pressure  drop  characteristics  in  tree-like  microchannel  nets  are  numerically  investigated  and 
compared  to  the  corresponding  characteristics  in  traditional  serpentine  flow  patterns,  by  solving  the  Navier-Stokes  and  energy  equation  for 
an  incompressible  fluid  with  constant  properties  in  three  dimensions.  A  constant  heat  flux  is  applied  to  the  walls  of  the  square  cross-sectional 
channels.  The  intrinsic  advantage  of  tree-like  nets  with  respect  to  both  heat  transfer  and  pressure  drop  is  demonstrated.  In  addition,  secondary 
flow  motions  initiated  at  bifurcations  and  their  important  role  on  thermal  mixing  are  identified  and  discussed.  Thermal  management  issues 
in  polymer  electrolyte  fuel  cells  are  addressed  and  in  this  context,  the  future  employment  of  tree  nets  is  recommended. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Effective  thermal  and  water  management  in  polymer  elec¬ 
trolyte  fuel  cells  is  key  to  ensure  high  cell  performance  and 
efficiency.  The  irreversibility  of  electrochemical  reactions 
and  joule  heating  are  the  most  important  factors  causing 
heat  generation  inside  PEM  fuel  cells.  The  temperature  dis¬ 
tribution  in  the  cell  has  a  strong  impact  on  the  cell  perfor¬ 
mance  since  it  influences  the  water  distribution  by  means  of 
condensation  and  affects  the  multicomponent  gas  diffusion 
transport  characteristics  through  thermocapillary  forces  and 
thermal  buoyancy.  Also,  the  kinetics  of  electrochemical  re¬ 
actions  directly  depends  on  temperature.  Excessive  global 
or  local  cell  temperature  due  to  insufficient  or  non-effective 
cell  cooling  may  cause  membrane  dehydration,  shrinking  or 
even  rupture.  Membrane  hydration  is  important  to  ensure 
high  proton  conductivity  and  thus  cell  performance.  Hence, 
thermal  and  water  management  issues  are  strongly  coupled 
and  they  have  a  direct  impact  on  cell  performance.  Ther¬ 
mal  management  includes  the  removal  of  the  generated  heat 
from  inside  the  cell  to  the  outside.  Further,  a  temporally  and 
spatially  uniform  temperature  distribution  must  be  provided, 
hot  spots  need  to  be  avoided,  and  pumping  power  required 
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for  the  coolant  circulation  has  to  be  minimized  in  order  to 
ensure  high  overall  cell  efficiency.  Therefore,  pressure  drop 
must  be  minimized  while  maximizing  the  heat  transfer  ca¬ 
pability  at  the  same  time. 

For  thermal  management  in  PEM  fuel  cell  stacks,  usu¬ 
ally  serpentine  or  meander  cooling  patterns  are  used,  hav¬ 
ing  only  one  or  a  few  inlets  and  outlets,  for  liquid  or  gas 
coolant,  depending  on  the  system  design  and  the  cell  size. 
For  overall  stack  functionality  reasons,  the  number  of  fluid 
inlets  and  outlets  must  be  small,  making  parallel  cooling 
patterns  usually  inappropriate.  These  circumstances  call  for 
a  flow  geometry  with  minimum  flow  resistance  between  a 
volume  and  one  point,  i.e.  between  the  inside  of  the  cell  and 
the  inlet  or  outlet.  In  this  paper,  the  use  of  deterministic  tree 
networks  of  Bejan  [1,2],  also  termed  as  fractal  networks  of 
branching  tubes  [3],  is  suggested  for  this  purpose.  Bejan  and 
Errera  [4]  discussed  a  deterministic  tree  network  for  fluid 
flow  with  minimum  flow  resistance  between  a  volume  and 
a  point,  subjected  to  two  constraints:  fixed  total  volume  and 
fixed  channel  volume.  Bejan  [5]  showed  that  the  total  heat 
current  convected  by  a  double  tree  is  proportional  to  the  to¬ 
tal  volume  raised  to  power  3/4.  Chen  and  Cheng  [6],  follow¬ 
ing  the  principles  put  forth  by  Bejan  [1,2,5],  presented  an 
analytical  study  on  the  heat  transfer  and  pressure  drop  char¬ 
acteristics  in  fractal  tree-like  microchannel  nets  in  which 
a  comparison  is  made  between  fractal  nets  and  traditional 
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Nomenclature 

a 

distance  (m) 

Br 

Brinkman  number  (-) 

cp 

specific  heat  capacity  (Jkg-1  K-1) 

d 

channel  diameter  (m) 

Ec 

Eckert  number  (-) 

f 

Moody  (or  Darcy)  friction  factor  (-) 

h 

heat  transfer  coefficient  (W  m-2  K-1) 

k 

branching  level  (-) 

L 

channel  length  (m) 

m 

total  number  of  branching  levels  (-) 

m 

mass  flow  rate  (kgs-1) 

n 

parameter  (-) 

Nu 

Nusslet  number  (-) 

P 

pressure  (Pa) 

Pr 

Prandtl  number  (-) 

q" 

wall  heat  flux  (W  m-2) 

r 

length  ratio  (-) 

Re 

Reynolds  number  (-) 

S 

surface  area  (m2) 

t 

time  (s) 

T 

temperature  (K) 

V 

velocity  (ms-1) 

V,  y,  z 

coordinates  (m) 

Greek  letters 

p 

channel  diameter  ratio  (-) 

r 

channel  length  ratio  (-) 

X 

thermal  conductivity  (W  m-1  K-1) 

ti 

dynamic  viscosity  (Pa  s) 

coordinates  (m) 

P 

density  (kgm-3) 

Subscripts 

Uj 

indices  (Einstein  summation  convention) 

in 

inlet 

k 

branching  level 

m 

mean 

max 

maximum 

min 

minimum 

ref 

reference 

w 

wall 

& 

in  §£  direction  where  k  indicates  the 
branching  level 

0 

zeroth  branching  level 

Superscript 

dimensionless 

parallel  channels.  After  simplifying  assumptions  including 
neglecting  the  effect  of  bifurcation  on  pressure  drop  and  heat 
transfer  they  found  that  the  considered  fractal  nets  provide 
higher  heat  transfer  capability  and  cause  less  pressure  drop. 


In  the  course  of  electronics  miniaturization  during  the 
last  decade,  denser  component  packing  led  to  higher  en¬ 
ergy  dissipation  rates  and  necessitated  the  demand  for 
high-performance  micro  heat  exchangers  consisting  of  mi- 
crochannel  systems  being  based  on  gas  or  liquid  convective 
heat  transfer  in  single  or  two-phase  flow  regimes  [7-9]. 
Compared  to  large  diameter  channels,  the  large  surface  to 
volume  ratio  of  small  diameter  channels  provides  a  better 
heat  transfer  capability  and  therefore  makes  microchannel 
systems  attractive  for  dealing  with  high  heat  fluxes.  How¬ 
ever,  a  penalty  is  paid  due  to  the  higher  pressure  drop, 
caused  again  by  the  increased  surface  to  volume  ratio.  In  the 
literature,  microchannels  are  generally  defined  as  channels 
having  a  hydraulic  diameter  smaller  than  1  mm  and  larger 
than  1  pan  [7]. 


Bipolar  plate,  anode  side 


(a)  Bipolar  plate,  cathode  side 


Fuel  cell  stack 


Gas  diffusion  layer,  cathode  side 

Membrane 
Gas  diffusion  layer,  anode  side 


Bipolar  plate,  anode  side 
Separation  layer 
Bipolar  plate,  cathode  side 


/ 


/ 

Cooling  layer  P 


Bipolar  plate,  anode  side 


Fig.  1.  Thermal  management  in  polymer  electrolyte  fuel  cell  stacks,  (a) 
Traditional  serpentine  or  meander  cooling  patterns  without  a  separation 
layer  between  the  bipolar  plates,  (b)  Superimposed  tree  nets.  One  net  is 
imposed  on  the  backside  of  the  anode  bipolar  plate  and  another  net  is 
imposed  on  the  backside  of  the  cathode  bipolar  plate.  The  two  nets  are 
separated  by  a  separation  layer.  Vertical  holes  in  this  layer  connect  the 
ends  of  the  highest  branching  levels  to  each  other. 
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A  substantial  research  effort  has  been  focused  on  the  ex¬ 
perimental  [10-17,19,21]  and  numerical  [18-24]  investiga¬ 
tion  of  single  phase  flow  in  such  microchannels.  Most  of 
this  work  described  in  literature  is  concerned  with  two  pa¬ 
rameters,  the  friction  factor  and  the  heat  transfer  coefficient, 
regarding  channels  with  diameters  between  1  jxm  and  1  mm. 
Experimental  results  show  agreement  but  also  substantial  de¬ 
viation  from  classical  theory  developed  for  larger  channels. 
Possible  reasons  for  the  measured  deviations  are  discussed 
and  investigated  in  the  literature  [25-33]:  measurement  in¬ 
accuracy  at  the  micro-scales,  high  relative  surface  roughness 
at  the  channel  walls,  the  electric  double-layer  (EDL)  effect, 
viscous  heating,  and  slip-flow  for  gases  at  high  Knudsen 
numbers.  Due  to  the  large  surface  to  volume  ratio,  surface 
related  phenomena  become  important  or  rather  dominant  at 
the  micro-scales  and  they  might  not  be  neglected,  as  at  the 
larger  scales.  However,  the  literature  is  not  conclusive  with 
respect  to  the  dependence  of  the  above  mentioned  two  pa¬ 
rameters  on  the  hydraulic  diameter  in  the  submicron- scales. 
A  key  property  of  a  tree  net  is  that  there  is  not  an  intrinsic 
length-scale,  defined  in  an  ad  hoc  manner.  Such  scales  are 
reduced  from  one  branching  level  to  the  next  higher  one, 
such  that  the  relative  importance  of  different  phenomena  de¬ 
pends  on  the  branching  level. 

In  this  paper,  a  numerical  analysis  of  the  fundamental 
heat  transfer  and  pressure  drop  characteristics  in  tree-like 
microchannel  nets  as  shown  in  Figs.  1  and  2  is  presented. 


The  future  application  of  tree  nets  for  thermal  management 
in  PEM  fuel  cells  is  suggested  and  discussed,  among  other 
application  areas.  In  the  analytical  work  by  Chen  and  Cheng 
[6],  the  effect  of  bifurcation  on  heat  transfer  and  pressure 
drop  is  neglected  and  it  is  assumed  for  simplicity,  that  the 
laminar  flow  is  fully  developed  both  thermally  and  hydro- 
dynamically.  Such  assumptions  are  appropriate  if  for  the 
flow  channels,  the  length  to  diameter  ratio  is  very  high  and 
the  Reynolds  number  is  rather  small,  such  that  the  hydro- 
dynamic  and  thermal  development  lengths  can  be  neglected 
compared  to  the  channel  lengths.  The  present  computational 
study  is  not  based  on  such  assumptions  and  it  takes  into  ac¬ 
count  important  hydrodynamic  and  thermal  effects  from  bi¬ 
furcations,  such  as  secondary  flow  motions  and  convective 
thermal  mixing.  The  Navier-Stokes  equations  and  the  en¬ 
ergy  equation  are  numerically  solved  in  three  dimensions  in 
conjunction  with  a  constant  heat  flux  boundary  condition  at 
the  channel  walls.  The  characteristic  variations  of  the  Nus- 
slet  number  and  the  friction  factor  along  the  tree  net  are  dis¬ 
cussed  and  compared  to  the  ones  of  a  serpentine  flow  pattern 
having  the  same  heat  transfer  surface  area. 

2.  Mathematical  model 

In  the  following  mathematical  model,  an  incompressible 
fluid  is  assumed  with  constant  viscosity  /x,  constant  heat 


Fig.  2.  Tree  net  with  six  branching  levels  m  =  6,  /3  =  <4+l/<4  =  2-1/3,  y  =  Lk+i/Lk  =  2-1//2,  Lo/do  =  10,  r  =  0.95,  and  A dk  =  a (a) 
Three-dimensional  view,  (b)  Top  view,  (c)  Schematic  of  the  geometric  structure  of  the  tree  net.  (d)  Computational  domain  with  inlet  and  outlet  mass 
flow  rates. 
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capacity  cp,  and  constant  thermal  conductivity  A.  Using  the 
Einstein  summation  convention,  continuity  reads 


9u/ 

dxj 


momentum  conservation  is  written  as 


3 

dt 


( PVi )  + 


a 

dxj 


( PVjVi ) 


dp  d 
dxi  dxj 


i=  1,2,3, 


and  the  energy  equation  including  the  viscous  dissipation 
term  is  given  by 


d  d 

-rXP1)  +  7, —  (PvjT) 
dt  dxj 

X  dT  \  p  /  3 Vi  di)j 
cp  dxj  )  cp  \  dxj  dxi 

For  a  non-dimensional  representation  of  the  problem,  the 
hydraulic  channel  diameter  do  of  the  zeroth  branching  level 
is  used  as  a  length-scale,  the  inlet  velocity  V[n  is  used  as 
a  velocity  scale,  the  inlet  temperature  Tm  and  a  maximum 
temperature  Tmax  are  used  for  temperature  scaling,  and  pres¬ 
sure  is  scaled  through  a  reference  pressure  pre f: 


3  Vi 


dx 


J 


(3) 


~ _  ^in 

~dk’ 


Vi 

Vi  =  - , 

Vm 


T-Tin 
Tmax  Tjn 


P  Pref 


The  maximum  temperature  is  defined  as  Tmax  =  q"  S/ 
( riik=ocp )  +  7in,  where  q"  is  the  constant  wall  heat  flux  and 
S  represents  the  corresponding  heated  surface  area.  Based 
on  these  scalings  and  the  definition  of  non-dimensional 
numbers, 


Re  = 


pvmdo 


Pr  = 


p 


c  pP 
~X~ 


Br  = 


X  ( Tmax  Tin) 


Ec  = 


Br 

Pr' 


(5) 


the  conservation  equations  (1)— (3)  can  be  written  in 
non-dimensional  form  as 


dv 


J 


dx 


=  0, 


j 


(6) 


3  „  3  _  dp  13 ^  Vi 

—  Vi~\ — —(vjVi)  = — —  H - - — — ,  f  =  1,2,3,  (7) 

ut  OX  j  oXf  uX  j  ox 


and 


3  ~  3 

—  T  H — — 
dt  dx 


CvjT) 


J 


1  3  ZT 


Re  Pr  dxjdxj 


+ 


Br 


3 Vi  3 vj\  dv 


Re  Pr  \  dx 


—  + 


J 


dx;  7  dx 


(8) 


j 


3.  Numerical  solution 

3.1.  Models  and  boundary  conditions 

Based  on  the  above  mathematical  formulation,  the  fluid 
mechanics  and  heat  transfer  in  a  tree-like  net  consisting  of 
channels  with  square  cross-sections,  as  shown  in  Fig.  2,  is 
computed.  The  fractal  geometry  considered  here  is  fixed  if 
the  following  independent  parameters  are  given:  the  number 
of  branching  levels  m,  the  ratio  of  the  channel  diameters 
/ 3  =  dk+i/dk  between  two  consecutive  branching  levels, 
the  corresponding  channel  length  ratio  y  =  Lk+i/Lk,  the 
length  to  diameter  ratio  Lo/do  at  the  zeroth  branching  level, 
the  ratio  r  of  the  straight  channel  length  and  the  entire  chan¬ 
nel  length  Lk  including  the  length  of  the  converging  duct, 
and  the  length  To  or  the  diameter  do  at  the  zeroth  branching 
level.  In  this  study,  a  fractal  geometry  with  m  =  6,  /3  = 
2-1/3,  y  =  2-1/2,  Lo/do  =  10,  r  =  0.95,  and  do  =  1mm 
is  investigated.  Optimization  of  a  single  plane  construct 
consisting  of  a  T-shaped  junction  with  respect  to  minimum 
flow  resistance  by  fixing  the  total  tube  volume  yields  the 
optimal  diameter  ratio  /3  =  dk+\/dk  =  2-1/3,  known  as 
Murray’s  law  [34,35],  which  is  independent  of  the  length 
ratio  y  =  Lk+i/Lk  and  geometry.  If  the  space  allocated 
to  the  construct  is  further  fixed,  i.e.  2Lk+\Lk  =  const.,  a 
second  minimization  of  the  flow  resistance  [5]  yields  the 
optimal  length  ratio  y  =  Lk+i/Lk  =  2-1/3.  However,  this 
length  ratio  creates  the  practical  problem  of  channel  over¬ 
lapping  for  a  right-angled  plane  construct  with  m  >  4.  To 
circumvent  this,  y  =  2_1//2  is  chosen  in  this  study.  Such  a 
tree  net  has  one  fluid  inlet  at  the  zeroth  branching  level  and 
a  number  of  2m  fluid  outlets  at  the  highest  branching  level, 
where  m  represents  the  total  number  of  branching  levels. 
A  traditional  heat  removal  concept  for  PEM  fuel  cells  is 
shown  in  Fig.  1(a)  in  which  a  serpentine  or  meander  flow 
pattern  is  manufactured  into  one  of  the  bipolar  plates. 

An  alternative  concept  utilizing  the  tree  network  of  chan¬ 
nels  requires  that  on  top  of  the  first  net,  a  second  net  of  iden¬ 
tical  geometry  is  superimposed  and  the  ends  of  the  highest 
branching  levels  of  each  net  are  then  vertically  connected  in 
order  to  build  up  a  closed  fractal  tree,  having  one  fluid  inlet 
at  the  zeroth  branching  level  of  the  first  net  and  one  fluid  out¬ 
let  at  the  zeroth  branching  level  of  the  second  net.  A  novel 
concept  for  thermal  management  in  PEM  fuel  cell  stacks 
which  is  based  on  this  principle  is  presented  in  Fig.  1(b): 
A  tree-like  flow  pattern  (first  tree  net)  is  imposed  on  the 
backside  of  the  anode  bipolar  plate  and  another  net  (second 
tree  net)  is  imposed  on  the  backside  of  the  adjacent  cathode 
bipolar  plate,  such  that  each  bipolar  plate  contains  a  flow 
pattern  for  reactant  and  product  gas  on  the  one  side  and  a 
tree-like  flow  pattern  for  thermal  management  on  the  other 
side.  Adjacent  bipolar  plates  are  separated  by  a  separation 
layer  including  vertical  holes  which  connect  the  ends  of  the 
highest  branching  levels  of  each  net.  In  this  manner,  the  an¬ 
ode  side  is  cooled  with  colder  coolant  whereas  the  cathode 
side  is  cooled  with  preheated,  warmer  coolant  providing  a 
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higher  temperature  level  on  the  cathode  side  than  on  the  an¬ 
ode  side.  Higher  temperatures  on  the  cathode  side  enhance 
electrochemical  reaction  rates  and  increase  the  saturation 
pressure  of  water  vapor  leading  to  a  reduced  amount  of  con¬ 
densation  which  in  turn  reduces  mass  transport  limitations 
at  high  current  densities.  On  the  cathode  side,  the  large  elec¬ 
trode  overpotential  and  the  condensation  of  water  vapor  are 
main  performance-limiting  factors  which  are  addressed  with 
the  proposed  concept.  Half  or  part  of  the  tree  net  may  be 
manufactured  on  the  bipolar  plate  and  the  other  half  on  the 
separation  layer.  Alternatively,  the  entire  tree  net  may  be 
manufactured  on  the  bipolar  plate  only.  Such  manufactur¬ 
ing  issues,  although  mentioned  in  passing  to  indicate  pos¬ 
sible  solutions,  are  not  the  main  theme  of  this  work  which 
focuses  on  the  basic  thermofluidic  aspects  of  the  problem. 
Heat  conduction  taking  place  in  the  solid  material  between 
the  flow  channels  (clearly  also  tractable  numerically)  is  not 
considered  in  this  work  opting  for  the  simple  constant  heat 
flux  condition  on  the  wall,  and  only  the  first  net  is  mod¬ 
eled  according  to  Fig.  2(a)-(c),  since  the  main  focus  is  to 
demonstrate  the  power  of  the  tree  channel  concept  leaving 
further  quantification  tied  to  specific  applications  for  future 
parametric  studies. 

This  net  is  symmetric  with  respect  to  the  planes  y  =  0 
and  z  =  0.  Therefore,  for  the  computational  domain  drawn 
in  Fig.  2(d),  symmetry  boundary  conditions  are  applied  at 
y  =  0  and  z  =  0.  A  constant  velocity  profile  is  applied  to  the 
inlet  boundary  condition.  The  tree  net  is  cut  at  the  left  end 
of  the  second  and  the  consequent  branching  levels  and  con¬ 
stant  mass  flow  outlet  boundary  conditions  are  used  to  re¬ 
duce  the  number  of  degrees  of  freedom  of  the  computational 
model,  cf.  Fig.  2(d).  This  configuration  of  outlet  boundary 
conditions  implicitly  incorporates  the  assumption,  that  at  a 
bifurcation,  half  of  the  mass  flow  is  directed  to  the  left  and 
the  other  half  to  the  right.  By  modeling  an  entire  tree  net 
with  m  =  4  for  a  wide  range  of  inlet  Reynolds  numbers,  it 
was  found  that  this  is  a  very  good  assumption  indeed,  even 
if  for  the  flow  it  is  not  possible  to  hydrodynamically  fully 
develop  before  bifurcating.  At  the  highest  branching  levels, 
fixed  pressure  outlet  boundary  conditions,  pre f  =  0,  together 
with  zero  flux  thermal  boundary  conditions  are  applied.  For 
heat  transfer,  a  constant  fluid  inlet  temperature  profile  is  set 
and  at  the  walls  of  the  tree  net,  a  constant  heat  flux  q"  and 
non-slip  boundary  conditions  are  applied. 

The  serpentine  geometry  as  shown  in  Fig.  3  has  exactly 
the  same  heat  transfer  surface  area  as  the  complete  frac¬ 
tal  geometry  with  an  identical  wall  heat  flux  q" .  Its  square 
cross-sectional  channels  are  of  constant  diameter  do  which 
equals  to  the  largest  channel  diameter  in  the  tree  net.  The 
entire  serpentine  flow  pattern  covers  the  same  rectangular 
area  as  the  tree  net. 

3.2.  Solution  method,  grid  independency 

The  set  of  partial  differential  equations  (6)-(8)  with  the 
boundary  conditions  described  above  is  numerically  solved 


Fig.  3.  Serpentine  flow  pattern  with  the  same  heat  transfer  surface  area 
as  the  tree  net  and  covering  the  same  rectangular  area. 

with  a  finite- volume  central  scheme  on  a  structured  grid  us¬ 
ing  the  licensed  solver  CFD-ACE  distributed  by  CFD  Re¬ 
search  Corporation  of  Huntsville,  USA.  Algebraic  multigrid 
solver  is  used  in  conjunction  with  a  convergence  criterion 
of  1  x  10-10. 

The  computational  grid  of  the  three-dimensional  fractal 
geometry  consists  of  1,816,800  cells.  By  increasing  the  cell 
number  to  2,721,705,  practically  indistinguishable  tempera¬ 
ture  contours  and  velocity  contours  of  all  three  components 
are  obtained  when  comparing  to  the  results  obtained  from 
the  coarser  grid.  By  decreasing  the  cell  number  to  790,425, 
small  discrepancies  are  obtained  for  temperature  and  veloc¬ 
ity  contours.  The  grid  resolution  of  the  serpentine  flow  ge¬ 
ometry  is  consistent  with  the  one  of  the  fractal  geometry 
and  it  contains  1,063,200  cells. 


4.  Validation 


The  numerical  results  are  validated  with  analytical  and  nu¬ 
merical  solutions  from  literature.  For  incompressible,  lami¬ 
nar  and  fully  developed  steady  flow  of  constant  viscosity  in 
a  square  duct,  an  exact  analytical  solution  [18]  can  be  ob¬ 
tained  for  the  mass  and  momentum  equations  (6)  and  (7). 
Within  this  framework,  the  exact  Moody  (or  Darcy)  friction 
factor/ [18]  is  given  as 


f  Re  =  24 


-i  -l 


oo 


1  —  19271  5  n  5  tanh  ^ ^ 


n  =  1,3, 


=  56.908. 


Iqbal  et  al.  [36]  analyzed  the  H2  boundary  condition  [18] 
of  rectangular  ducts  with  all  four  walls  heated  for  a  ther¬ 
mally  and  hydrodynamically  fully  developed  flow  with  con¬ 
stant  fluid  properties  and  neglecting  viscous  dissipation.  For 
a  square  duct,  they  obtained  Nu  =  3.091  by  means  of  a 
variational  approach.  Chandrupatla  and  Sastri  [37]  obtained 
A Ju  =  3.095  which  is  in  agreement.  Their  finite  difference 
solutions  were  computed  by  using  the  iterative  extrapolated 
Liebmann  method. 
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For  Rek=0  =  20,  the  flow  becomes  fully  developed  both 
thermally  and  hydrodynamically  in  the  first  channel  of  the 
serpentine  flow  pattern.  The  present  numerical  solution  of 
Eqs.  (6)-(8),  neglecting  viscous  dissipation  in  the  energy 
equation  like  in  the  references  mentioned  earlier  to  facilitate 
the  comparison,  yields  /  Re  =  56.976  (+0.118%  compared 
to  Ref.  [18])  and  Nu  =  3.096  (+0.162  and  +0.0323%  com¬ 
pared  to  Refs.  [36,37],  respectively),  which  is  in  very  good 
agreement  with  the  previously  mentioned  results. 


5.  Nusslet  number  and  friction  factor 


In  order  to  investigate  the  heat  transfer  and  pressure  drop 
characteristics  of  the  tree  net  and  the  serpentine  flow  pat¬ 
tern,  the  local  Nusslet  number  and  local  friction  factor  are 
determined  at  each  discrete  channel  cross-section  along  the 
flow  direction.  The  local  Nusslet  number  then  depends  on 
the  position  ^  as  indicated  in  Figs.  2(c)  and  3,  i.e. 


A Tuk(£k)  = 


hk(^k)dk 


(10) 


where  the  subscript  k  indicates  the  branching  level  in  the 
case  of  the  tree  net.  With  the  constant  wall  heat  flux  q ",  the 
local  heat  transfer  coefficient  reads 


M&)  = 


7w,m  ,ki^k) 


Tm,k(^k) 


(ID 


where  the  mean  fluid  temperature  is  given  as 


Tm,k(Hk) 

m,  ft)  dm  da 


So  ++  v^k’k ^k 'Vkitk)  dm  d a 


(12) 


and  the  peripherally  averaged  wall  temperature  reads 


7\v,m ,ki^k)  —  i  /nt  [  ,k(Hk>  0k,mm->  ^k)  dfp 

dk/^Jo 


ft, 


max 


+ 


+ 


|  fkk,  max 

“7  TZ  I  Tw,k(^k,  0k  .  max  •>  a)  d  a 
dk/2Jo 

^  f  Vk,  max 

I  Tw,k(^ki  Oki  Kk  ,max  )dm-  (13) 

Ok  J  m 


Vk,  min 


Herein,  v^k  is  the  velocity  component  in  ^-direction  at 
the  kth  branching  level.  The  Reynolds  number  at  the  kth 
branching  level, 

PV^k,  vafdk 


Rek  = 


(14) 


with  the  mean  velocity 


v%k, — 


fo  fnflTf  vHkM^k)  0k)  &)  dllk  d^k 


Vk,  min 

Ik, 


/o  C,n“d%dft 

°  w  °  'Ik, min 


=  2~kB~2kVin, 


(15) 


and  the  definition  of  the  friction  factor, 


/*(&)  = 


dpm,k  2  dp 

^k  f 


(16) 


are  used  to  compute  the  local  fpRepi^k)  which  depends  on 
the  position  along  the  flow  direction.  The  cross-sectional 
average  of  the  static  pressure 


Pm,k(^k)  — 


fo  Pk  ^0k  dfk 


Vk,  min 

Ik, 


foiled  md& 

J  u  J  'Ik, mm 


(17) 


is  used  to  determine  the  pressure  gradient  in  Eq.  (16). 


6.  Numerical  results  and  discussion 

The  equations  governing  fluid  flow  and  heat  transfer, 
Eqs.  (6)-(8),  are  numerically  solved  with  respect  to  the 
tree  and  serpentine  channel  patterns.  For  each  configura¬ 
tion,  two  cases  are  considered:  The  first  case  is  based  on 
Rek= o  =  20  and  Pr  =  2.39  whereas  the  second  case  is 
based  on  Rek= o  =  200  and  Pr  =  2.39.  This  Prandtl  number 
represents  the  properties  of  liquid  water  at  1  bar  and  75  °C 
[38].  The  operating  temperature  range  for  PEM  fuel  cells 
is  usually  between  65  and  100  °C.  For  a  channel  diameter 
of  do  =  1  mm  at  the  zeroth  branching  level,  these  Reynolds 
numbers  then  imply  a  mean  liquid  water  inlet  velocity  of 
0.774  and  7.74  cm/s,  respectively.  In  this  section,  laminar 
mixing,  pressure  drop,  and  heat  transfer  characteristics  in 
the  fractal  and  serpentine  flow  pattern  are  discussed  and 
compared  to  each  other.  Dimensional  results  are  based  on 
the  parameters  given  in  Table  1.  Characteristics  are  plotted 
for  the  right-sided  branches,  cf.  Fig.  2(d). 

6.1.  Laminar  mixing 

As  the  flow  passes  through  bifurcations,  secondary  flow 
motions  are  initiated  as  shown  in  Fig.  4(a)  and  (b).  The 
higher  the  flow  Reynolds  number,  the  larger  their  relative 
magnitude  is.  Transverse  vortices  with  their  axis  transverse 
to  the  flow  in  the  streamwise  direction  generate  recirculation 


Table  1 


Properties  of  liquid  water  at  1  bar  and  75  °C  [38].  Boundary  conditions 
and  geometry  parameters 


Density,  p  (kgm-3) 

975 

Dynamic  viscosity,  pi  (Pas) 

377  x  10-6 

Thermal  conductivity,  X  (Wm-1  K-1) 

0.664 

Specific  heat  capacity,  cp  (Jkg-1  K-1) 

4203 

Inlet  velocity  v$  for  Rek= o  =  20  (cms-1) 

0.774 

Inlet  velocity  v$  for  Rek= o  =  200  (cms-1) 

7.74 

Wall  heat  flux  q"  for  Rek= o  =  20  (Wm-2  K-1) 

776.8 

Wall  heat  flux  q"  for  Rek= 0  =  200  (Wm“2  K"1) 

7768 

Channel  diameter  do  at  the  zeroth  branching 
level  (mm) 

1.0 

Surface  area,  S  (m2) 

4.083  x  10“4 
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Longitudinal 

Vortices 


Fig.  4.  Laminar  mixing  at  the  first  bifurcation  in  the  tree  net  at 
Rek= o  =  200.  (a)  Streamlines  in  top  view:  transverse  vortices  with  their 
axis  transverse  to  the  flow  causing  recirculation  and  longitudinal  vor¬ 
tices  with  their  axis  along  the  flow  in  the  streamwise  direction,  (b) 
Three-dimensional  view,  (c)  Temperature  (K)  contours  showing  hot  spots 
in  the  recirculation  area  at  z  =  0. 


near  the  inner  corners  of  the  bifurcations,  where  the  fluid 
velocity  is  much  slower  than  at  the  outer  wall,  leading  to 
hot  spots  in  this  area,  cf.  Fig.  4(c).  Further,  two  longitudinal 
vortices  of  opposite  orientation  with  their  axis  along  the  flow 
in  the  streamwise  direction  are  generated,  as  seen  in  Fig.  5. 
Thereby,  cold  fluid  is  moved  from  the  outer  to  the  inner  wall 
and  finally  to  the  channel  core  as  the  strength  of  the  vortices 
decays  with  increasing  distance  from  the  bifurcation,  cf. 
Fig.  5.  With  increasing  branching  level  k ,  the  flow  Reynolds 
number  becomes  smaller  and  thereby  the  initial  strength  of 
longitudinal  vortices  is  reduced  such  that  cold  fluid  is  not 
anymore  subject  to  establish  a  closed  ring  with  warmer  fluid 
inside  this  ring. 

For  the  tree  net,  the  distance  A dk  as  indicated  in  Fig.  2(c) 
has  been  varied  to  investigate  the  dependence  of  the  vor¬ 
tex  structures  on  different  geometric  configurations  of  the 
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Fig.  5.  (a)  Temperature  (K)  contours  at  four  consecutive  cross-sections 
in  the  first,  second,  and  third  branching  level  of  the  tree  net  at 
Rek= o  =  200.  (b)  Vorticity  (s_1)  and  pressure  (Pa)  contours  at  two  con¬ 
secutive  cross-sections  in  the  first  branching  level  k  =  1  of  the  tree  net 
at  Rek= o  =  200. 
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converging  passages.  Identical  longitudinal  vortex  struc¬ 
tures  of  varying  strength  are  found  for  different  values  of 
A dk.  Herein,  the  following  parameter  set  has  been  investi¬ 
gated  for  the  higher  Reynolds  number  level  (Rek= o  =  200) : 
A dk  €  (-2 ak,  -ak,0,ak,2ak)  with  ak  =  (dk  -  dk+ 1)/2 
where  the  geometry  with  A dk  =  ak  represents  the  case 
being  mainly  discussed  in  this  study.  The  strength  of  the 
longitudinal  vortices  increases  for  larger  values  of  A dk,  in¬ 
creasing  thereby  the  peak  local  Nusslet  number  Nuk=  i?max 
according  to  10.2,  10.9,  12.7,  15.1,  and  17.0.  However,  the 
pressure  drop  in  the  entire  tree  net  increases  at  the  same  time 
according  to  50.2,  52.1,  56.1,  65.0,  and  87.3  Pa,  respectively. 
Therefore,  the  effect  of  thermal  mixing  can  be  intensified 
with  increasing  values  of  A dk  at  the  cost  of  additional  pres¬ 
sure  drop  and  vice  versa.  The  channel  cross-sectional  area 
at  the  interface  of  two  consecutive  branching  levels  k  and 
k  +  1,  i.e.  ( dk  —  2Adk)dk+i,  is  reduced  with  increasing  A dk. 


Therefore,  the  flow  is  locally  accelerated  before  entering 
the  higher  branching  level  which  causes  higher  strengths 
for  longitudinal  vortices  and  therefore  enhanced  thermal 
mixing.  More  detailed  optimization  of  the  bifurcation  ge¬ 
ometry  at  the  different  branching  levels  and  inlet  Reynolds 
number  with  respect  to  both  aspects  is  left  to  algorithmic 
optimization  studies  considering  also  complexly  curved 
walls.  Fluid  mixing  could  also  be  intensified  by  different 
passive  methods,  e.g.  by  placing  ridges  on  the  floor  of  the 
channels  at  an  oblique  angle  [39]  or  by  inserting  C-shaped 
repeating  units  in  the  streamwise  direction  [40]. 

6.2.  Pressure  drop 

In  the  tree  net,  the  flow  is  subject  to  develop  hydrody- 
namically  and  thermally  in  the  entrance  region  of  the  zeroth 
branching  level  due  to  the  constant  velocity  and  temperature 


(d)  Position 


Fig.  6.  Friction  factor,  (a)  Tree  net  with  Rek= o  =  20.  (b)  Serpentine  pattern  with  Rek= o  =  20.  (c)  Tree  net  with  Rek= o  =  200.  (d)  Serpentine  pattern  with 
Rek= o  =  200. 
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Fig.  7.  Pressure  drop  in  the  straight  channels  between  ^=o,min  and  £&=o,max  (x).  Pressure  drop  across  bifurcations  between  £^max  and  ^+i,min  (+)•  (a) 
Tree  net  with  Rek=o  =  20.  (b)  Tree  net  with  Rek= o  =  200. 


profiles  at  the  inlet.  At  bifurcations,  the  flow  is  disturbed 
and  secondary  flow  motions  are  initiated  which  then  decay 
along  the  next  following  straight  channel  portions,  such  that 
the  flow  tends  to  develop  again  before  reaching  the  next  bi¬ 
furcation.  In  the  serpentine  channels,  the  flow  is  subject  to 
develop  in  the  straight  channels  whereas  it  is  disturbed  again 
at  each  turn. 

For  Rek=0  =  20,  the  hydrodynamic  development  lengths 
are  small  compared  to  the  length  of  the  corresponding 
straight  channels.  This  is  seen  in  Fig.  6(a)  and  (b)  from 
the  course  of  the  friction  factor  which  evens  out  at  a 
certain  level  when  the  flow  becomes  fully  developed.  In 
the  tree  net,  the  pressure  drop  across  the  first  bifurcation 
(0.145  Pa)  is  16.2%  of  the  pressure  drop  between  §&=o,min 
and  §£=o,max  in  the  straight  channel  of  the  zeroth  branching 
level  (0.898  Pa),  cf.  Fig.  7(a).  Herein,  the  cross-sectional 
average  of  the  static  pressure  at  ^max  is  subtracted  from 
the  corresponding  value  at  ^+i,min  to  evaluate  the  pressure 
drop  across  a  bifurcation.  The  pressure  drop  across  the 
sixth  bifurcation  (0.1 14  Pa)  is  28.3%  of  the  pressure  drop  in 
the  straight  channel  of  the  fifth  branching  level  (0.404  Pa). 
Consequently,  at  the  higher  branching  levels,  the  contri¬ 
bution  of  the  pressure  drop  obtained  across  bifurcations 
becomes  more  important  compared  to  the  pressure  drop  in 
the  straight  channels,  as  seen  from  Fig.  7(a).  This  fact  is  due 
to  the  growing  ratio  between  the  channel  diameter  and  the 
channel  length  when  going  to  higher  branching  levels,  i.e. 
dk+i/Lk+i  =  0 6/y)dk/Lk  =  l.\2dk/Lk.  Therefore,  analyt¬ 
ical  approaches  [6]  where  the  laminar  flow  is  assumed  to  be 
hydrodynamically  fully  developed  everywhere  and  where 
the  effect  of  bifurcation  on  pressure  drop  is  neglected,  are 
only  adequate  if  the  hydrodynamic  development  lengths  are 
very  small  compared  to  the  length  of  the  straight  channels, 
i.e.  for  low  Reynolds  numbers,  high  L^/dk  ratios  and  for 
the  lower  branching  levels.  For  the  entire  tree  net,  the  re¬ 
sulting  overall  pressure  drop  (4.88  Pa)  is  only  56.9%  of  the 


pressure  drop  obtained  from  the  corresponding  serpentine 
flow  pattern  (8.57  Pa)  having  exactly  the  same  surface  area 
S.  As  a  further  comparison,  the  pressure  drop  in  a  single 
square  cross-sectional  straight  channel  with  diameter  do  and 
surface  area  S  can  be  obtained  analytically  from  Eqs.  (9) 
and  (16)  by  assuming  fully  developed  flow.  A  value  of 
8.48  Pa  is  found  for  such  a  straight  channel  flow  which  is 
only  slightly  below  the  pressure  drop  in  the  serpentine  flow 
pattern  but  at  the  same  time  almost  double  the  pressure 
drop  in  the  tree  net. 

For  Rek= o  =  200,  the  flow  does  not  fully  develop  any 
more  at  the  lower  branching  levels  k  =  0,  1,  2  of  the  tree 
net,  as  seen  from  Fig.  6(c).  In  the  serpentine  pattern,  the  flow 
develops  in  the  long  straight  channels,  k  =  0,  2,  4,  such  that 
the  friction  factor  variations  in  the  shorter  channels  k  =  1 
and  k  —  3  are  similar  to  each  other,  cf.  Fig.  6(d).  The  pres¬ 
sure  drop  across  the  first  bifurcation  in  the  tree  net  (4.03  Pa) 
is  31.5%  of  the  pressure  drop  in  the  straight  channel  of  the 
zeroth  branching  level  (12.8  Pa),  cf.  Fig.  7(b).  The  pressure 
drop  across  the  fifth  bifurcation  (1.29  Pa)  is  25.4%  of  the 
pressure  drop  in  the  straight  channel  of  the  fourth  branching 
level  (5.09  Pa).  Here,  the  relative  importance  of  the  bifur¬ 
cation  pressure  drop  does  not  monotonically  increase  with 
increasing  branching  levels,  since  the  flow  does  not  develop 
fully  at  the  lower  branching  levels  leading  to  a  higher  rela¬ 
tive  pressure  drop  from  bifurcations  at  these  branching  lev¬ 
els.  For  k  =  0  to  k  =  5,  these  ratios  are  31.5,  26.9,  24.8, 
24.5,  25.4,  and  27.1%,  respectively.  As  expected,  the  ratios 
increase  again  after  the  third  branching  level  where  the  flow 
tends  to  develop  more  approaching  fully  developed  state  be¬ 
fore  bifurcating.  The  overall  pressure  drop  resulting  from 
the  entire  tree  net  (65.0  Pa)  is  only  55.1%  of  the  pressure 
drop  obtained  from  the  corresponding  serpentine  flow  pat¬ 
tern  (11 8  Pa).  When  comparing  the  overall  pressure  drop  of 
the  tree  net  at  Re^~ o  =  20  (4.88  Pa)  to  the  one  at  Re^~ o  = 
200  (65.0  Pa),  an  increase  by  a  factor  of  13.3  is  identified. 
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For  fully  developed  flow  and  neglecting  effects  from  bifur¬ 
cations,  a  factor  of  10  would  be  expected.  When  considering 
again  a  single  straight  channel  with  fully  developed  flow, 
a  pressure  drop  of  84.8  Pa  is  found  analytically,  which  is 
substantially  below  the  pressure  drop  in  the  serpentine  flow 
pattern  and  at  the  same  time  also  substantially  above  the 
pressure  drop  in  the  tree  net. 

In  conclusion,  the  relative  effect  of  pressure  drop  from 
bifurcations  becomes  more  dominant  at  the  higher  branch¬ 
ing  levels.  In  spite  of  the  large  number  of  bifurcations  in 
the  tree  net  compared  to  a  much  lower  number  of  turns  in 
the  serpentine  flow  pattern,  the  tree  net  proves  significantly 
beneficial  in  terms  of  pressure  drop. 

6.3.  Heat  transfer 

In  the  tree  net,  the  flow  becomes  thermally  fully  devel¬ 
oped  at  each  branching  level  for  Rek= o  =  20.  This  is  seen 


in  Fig.  8(a)  from  the  variation  of  the  Nusslet  number  which 
evens  out  at  a  certain  level  when  the  flow  is  subject  to  de¬ 
velop.  The  same  applies  to  the  serpentine  flow  pattern,  cf. 
Fig.  8(b).  As  expected,  the  heat  transfer  coefficient  hk  in 
the  thermally  fully  developed  region  scales  with  hk  oc  l/<4 
which  implies  that  in  the  highest  level  of  the  tree  net  a 
more  than  four  times  higher  heat  transfer  coefficient  can  be 
obtained  than  in  the  serpentine  flow  pattern,  as  seen  from 
Fig.  9(a)  and  (b).  Therefore,  a  smaller  temperature  differ¬ 
ence  is  needed  to  transfer  the  same  heat  flux  from  the  chan¬ 
nel  walls  to  the  bulk  fluid  at  the  higher  branching  levels  of 
the  tree  net.  Due  to  the  constant  wall  heat  flux  boundary 
condition  used  in  this  study  and  since  the  channel  diame¬ 
ter  gets  smaller  when  going  to  higher  branching  levels,  the 
length-specific  heat  transfer  rate  to  the  fluid  decreases  from 
one  branching  level  to  the  next  higher  one,  i.e.  by  the  factor 
ft  =  0.794.  The  heat  transfer  rate  to  a  single  channel  of  the 
branching  level  k  +  1  is  y/3  =  0.561  times  smaller  than  the 


Fig.  8.  Nusslet  Number,  (a)  Tree  net  with  Rek= o  =  20.  (b)  Serpentine  pattern  with  Rek= o  =  20.  (c)  Tree  net  with  Rek= o  =  200.  (d)  Serpentine  pattern 
with  Rek= o  —  200. 
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Fig.  9.  Heat  transfer  coefficient,  (a)  Tree  net  with  Reu= o  =  20.  (b)  Serpentine  pattern  with  Re a=o  =  20.  (c)  Tree  net  with  Re^= o  =  200.  (d)  Serpentine 
pattern  with  Re^= o  —  200. 


heat  transfer  rate  to  a  single  channel  of  the  branching  level 
k.  However,  the  heat  transfer  rate  to  all  the  channels  of  the 
branching  level  k  +  1  is  2 y/3  =  1.12  times  larger  than  the 
heat  transfer  rate  to  all  the  channels  of  the  branching  level 
k.  Only  for  large  Lk/dk  ratios  and  low  Reynolds  numbers 
or  if  the  thermal  development  length  can  be  neglected  with 
respect  to  the  channel  length  L^,  the  assumption  of  constant 
Nusslet  number  to  compute  the  heat  transfer  capability  of 
the  tree  net  is  adequate. 

For  a  10  times  higher  inlet  Reynolds  number  to  the 
tree  net,  i.e.  Re^— o  =  200,  the  flow  does  not  thermally 
develop  any  more,  except  at  the  highest  branching  levels, 
as  seen  from  Fig.  8(c).  When  the  flow  passes  through  the 
bifurcation,  secondary  flow  motions  are  initiated  that  de¬ 
generate  the  thermal  boundary  layer.  This  effect  causes 
high-temperature  fluid,  which  is  captivated  near  the  walls 
in  the  fully  developed  region,  to  mix  with  low-temperature 
fluid  from  the  channel  core,  as  illustrated  in  Fig.  5.  This  con¬ 
vective  thermal  mixing  effect  provides  a  smaller  difference 


between  the  mean  wall  temperature  Tw?m ^  and  the  mean 
fluid  temperature  Tm ^  and,  due  to  the  constant  wall  heat 
flux,  leads  to  an  increase  of  the  local  Nusslet  number  Nuk 
in  the  straight  channels  right  after  bifurcations,  cf.  Fig. 
8(c).  The  effect  of  thermal  mixing  caused  by  secondary 
flow  motions  and  reinitiated  at  each  bifurcation  becomes 
dominant  at  this  inlet  Reynolds  number  level  and  the  Nus¬ 
slet  number  reaches  a  peak  value  of  Nuk=\  =  15.1  at 
a  distance  §£=i/§£=i,max  =  0.106  after  bifurcating.  This 
value  is  almost  five  times  larger  than  the  Nusslet  number 
of  a  thermally  and  hydrodynamically  fully  developed  flow. 
For  higher  branching  levels,  the  peak  Nusslet  number  de¬ 
creases  from  one  level  to  the  next  higher.  The  effectiveness 
of  thermal  mixing  initiated  at  bifurcations  depends  on  the 
Reynolds  number.  It  is  improved  for  increasing  Reynolds 
numbers  and  vice  versa.  Since  the  Reynolds  number  de¬ 
creases  from  one  branching  level  to  the  next  higher  one,  i.e. 
according  to  Re^x/Re^  =  l/(2/3)  =  0.630,  thermal  mixing 
becomes  less  effective  at  the  higher  branching  levels  which 
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manifests  itself  in  the  above-mentioned  decrease  of  the  peak 
Nusslet  number  from  one  level  to  the  next  higher  one.  The 
local  heat  transfer  coefficient  increases  from  one  branching 
level  to  the  next  higher  with  high  local  peaks  due  to  thermal 
mixing,  cf.  Fig.  9(c).  In  many  applications,  e.g.  fuel  cell 
cooling  and  electronics  cooling,  a  thermofluidic  design  for 
thermal  management  is  subject  to  the  constraint  of  a  max¬ 
imum  temperature  difference  within  the  considered  system. 
Therefore,  when  attempting  to  transfer  larger  heat  fluxes 
under  such  a  constraint,  higher  heat  transfer  coefficients  are 
required  since  they  allow  a  certain  heat  flux  to  be  transferred 
at  a  smaller  temperature  difference.  In  this  respect,  the  tree 
net  proves  significantly  beneficial,  as  shown  in  Fig.  9(a)-(d). 

The  higher  the  branching  level,  the  more  beneficial  it 
proves.  Hence,  when  considering  a  single  tree  net,  a  higher 
maximum  heat  rate  could  be  transferred  than  with  the  ser¬ 
pentine  flow  pattern,  if  a  certain  maximum  temperature 
difference  is  imposed  on  the  two  systems,  as  seen  from 


Fig.  10(a)— (d).  On  the  other  hand,  reduced  temperature  dif¬ 
ferences  needed  to  transfer  a  constant  wall  heat  flux  at  the 
higher  branching  levels  imply  a  reduced  maximum  temper¬ 
ature  difference  within  the  entire  system.  At  Re^— o  =  200, 
a  maximum  temperature  difference  between  rw?max  and 
Tm,min  of  12.0  K  is  obtained  from  the  tree  net,  compared  to 
the  corresponding  value  of  16.0  K  from  the  serpentine  flow 
pattern.  For  Rek= o  =  20  and  with  a  10  times  lower  wall 
heat  flux  than  in  the  higher  Reynolds  number  case,  these 
two  values  are  closer  together,  as  seen  from  Fig.  10(a)-(d). 
The  higher  the  wall  heat  flux  and  the  larger  the  total  number 
of  branching  levels,  the  more  advantageous  the  single  tree 
net  with  respect  to  small  temperature  differences.  Providing 
a  uniform  temperature  distribution  is  one  of  the  previously 
mentioned  important  requirements  for  effective  thermal 
management  in  polymer  electrolyte  fuel  cells. 

For  a  double  tree  net,  this  argumentation  no  longer  holds. 
The  same  maximum  temperature  difference  will  be  obtained 


Fig.  10.  Mean  fluid  temperature,  mean  wall  temperature,  and  maximum  wall  temperature,  (a)  Tree  net  with  Rek= o  =  20,  q"  =  776.8  Wm  2K  1 . 
(b)  Serpentine  pattern  with  Rek= o  =  20,  q "  =  776.8  Wm-2  K-1.  (c)  Tree  net  with  Rek= o  =  200,  q"  =  7768Wm-2K-1.  (d)  Serpentine  pattern  with 
Rek= o  =  200,  q"  =  7768  Wm^K'1. 
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from  a  double  tree  net  and  from  a  serpentine  flow  pattern 
of  identical  surface  area  and  with  identical  constant  wall 
heat  flux.  However,  it  is  still  thermodynamically  desirable 
to  transfer  a  certain  heat  flux  with  a  small  temperature  dif¬ 
ference  because  of  smaller  corresponding  entropy  genera¬ 
tion.  In  addition,  there  is  room  for  the  transfer  of  higher  heat 
fluxes  if  needed,  by  increasing  the  temperature  difference. 
As  an  example,  a  substantially  higher  overall  heat  transfer 
rate  could  be  obtained  from  the  single  or  double  tree  net  com¬ 
pared  to  the  serpentine  flow  pattern  of  the  same  surface  area, 
if  instead  of  the  constant  wall  heat  flux  a  constant  tempera¬ 
ture  difference,  rw?m^  —  Tm^  =  const.,  could  be  maintained 
along  the  entire  channel  system.  This  can  be  seen  from  the 
variation  of  the  mean  wall  temperature  and  the  mean  fluid 
temperature  in  Fig.  10(a)  and  (b):  In  the  fully  developed  re¬ 
gions  of  the  tree  net,  the  difference  between  the  mean  wall 
temperature  and  the  mean  fluid  temperature  decreases  from 
one  branching  level  to  the  next  higher  whereas  in  the  ser¬ 
pentine  pattern  this  difference  remains  constant.  Enforcing 
such  a  constant  temperature  difference  in  all  the  branching 
levels  of  the  tree  net  would  imply  larger  heat  fluxes  at  the 
higher  branching  levels  and  therefore  the  curves  indicating 
the  variation  of  the  mean  fluid  temperature  would  get  steeper 
from  one  branching  level  to  the  next  higher.  Hence,  under 
such  conditions,  a  higher  mean  fluid  temperature  would  be 
obtained  at  the  end  of  the  three  net  than  at  the  end  of  the 
serpentine  pattern.  Alternatively,  in  trying  to  get  the  same 
mean  fluid  temperature  at  the  end  of  both  configurations,  a 
tree  net  with  a  smaller  number  of  branching  levels  would  be 
necessary  which  would  further  reduce  the  amount  of  pres¬ 
sure  drop  in  the  tree  net. 

Finally,  it  is  worth  mentioning  that  the  effect  of  viscous 
dissipation  can  be  safely  neglected  in  the  cases  considered 
so  far,  as  seen  from  the  local  Brinkman  numbers,  which  are 
in  the  order  of  1  x  10-6.  This  was  also  numerically  verified 
by  simply  dropping  the  viscous  dissipation  term  out  of  the 
energy  equation  and  then  comparing  with  the  former  results. 

7.  Conclusions 

A  thorough  three-dimensional  forced  convection  model 
was  presented  for  two  geometrical  channel  configurations 
with  constant  wall  heat  flux:  a  tree-like  net  and  a  serpentine 
flow  pattern  having  the  same  surface  area.  The  tree  network 
configuration  represents  the  main  novel  concept  of  this  paper 
for  cooling  systems  in  applications  such  as  fuel  cells  and 
electronic  devices.  The  serpentine  configuration  is  typical 
in  the  cooling  of  PEM  fuel  cells  and  provides  a  basis  of 
comparison. 

It  is  found  that  the  tree  net  with  six  branching  levels  gen¬ 
erates  almost  half  the  pressure  drop  than  the  corresponding 
serpentine  flow  pattern  having  the  same  surface  area  and  in¬ 
let  Reynolds  number.  The  effect  of  pressure  drop  from  bifur¬ 
cations  and  turns  is  substantial  and  should  not  be  neglected 
for  simplicity  by  using  fully  developed  flow  friction  factors 


in  the  considered  cases.  Although  the  number  of  bifurcations 
in  the  tree  net  is  larger  than  the  number  of  turns  in  the  ser¬ 
pentine  flow  pattern,  from  which  one  would  expect  a  lower 
pressure  drop  from  the  serpentine  flow  pattern,  the  tree  net 
proves  to  be  markedly  beneficial  in  terms  of  pressure  drop. 

In  addition  to  hydrodynamic  advantages,  the  tree  net  fur¬ 
ther  provides  a  larger  heat  transfer  capability  than  the  ser¬ 
pentine  pattern.  Laminar  mixing  by  secondary  flow  motions 
initiated  at  bifurcations  causes  substantially  improved  local 
Nusslet  numbers.  Much  higher  heat  transfer  coefficients  are 
obtained  for  the  tree  net  which  implies  that  at  the  higher 
branch  levels,  a  certain  heat  flux  can  be  transferred  with  a 
smaller  temperature  difference.  A  much  higher  heat  transfer 
rate  could  be  obtained  for  the  tree  net  by  increasing  this  tem¬ 
perature  difference.  If  for  the  water  temperature  it  is  possible 
to  reach  the  boiling  temperature,  i.e.  if  the  cell  would  be  op¬ 
erated  above  this  temperature  level,  evaporation  can  occur, 
leading  to  different  regimes  of  two-phase  flow  in  the  cool¬ 
ing  channels.  This  case  was  not  investigated  in  the  present 
study  which  focused  on  the  single-phase  flow  regime. 

Due  to  their  intrinsic  advantage  with  respect  to  both  heat 
transfer  and  pressure  drop,  tree-like  nets  have  the  poten¬ 
tial  to  significantly  aid  the  thermal  management  in  polymer 
electrolyte  fuel  cells  at  different  length  scales  as  well  as  in 
other  thermal  systems. 
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